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Abstract 

The time-development of photoexcitations in the coupled chromophores exhibits specific dy- 
namics of electronic sites population and nuclear wavefunction. In many cases, the specifics of 
the site-population and wavefunction amplitude dynamics is determined by the initial state of the 
nuclear subsystem. We discuss the scenario of measuring the wavefunction of the system by means 
of nonlinear wavepacket interferometry that characterizes the dynamical entanglement formation 
of the vibronic quantum system in a consistent manner as opposed to the traditional population 
kinetics measurements. 
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I. INTRODUCTION 



These studies are motivated by the fundamental interest in properties of few-level elec- 
tronic system coupled to many-mode field, e.g. nuclear vibrations. The major goal is to 
check the possibility of controlling electronic population dynamics by varying the amount of 
vibrational excitation. Another goal is to characterize the state of vibronic quantum system 
on the amplitude level by means of ultrafast spectroscopy with account of the phase infor- 
mation^. The studies were also motivated by the potential application of this investigation 
to prospective systems, containing coupled chromophores e.g. natural and artificial light 
harvesting^ photographic imaging, and optical communication technologies^. 

The energy transfer pathway of electronic excitation of molecular systems has deserved 
a lot of attention last years^'^'^. In coupled chromophores system the dipole-dipole coupling 
J promotes the excitation from one chromophore to its neighbor sites^. The coherency 
between neighbor sites is usually destroyed due to the electronic-nuclear coupling character- 
ized by reorganization energy A. Various values of relation J/A embodies various regimes 
of dynamics^. Preparation of a nuclear mode in a specific state also affects the regime of 
system dynamics-'^'^. 

The progress in laser technology allows for direct measurements of the dynamical fea- 
tures of molecular systems by means of short-pulse spectroscopy. Among various suitable 
time-resolved spectroscopic techniques one could mention time-dependent fluorescent (also 
polarization resolved^^'^^'^^) measured by either photon counting or fluorescence upconver- 
sion technique^^'^^. An alternative option is pump-probe measurement of transient absorp- 
tion^^. The signal both methods is smeared out by inhomogeneity of molecular systems. 
The known technique that beats the inhomogeneity is photon echo. Three-pulse echo gives 
more information and allows to trace the vibrational dynamics in the population period 
between pulses. Homodyne (with the 4th pulse) detection is even more convenient because 
of one technical reason: Fluorescence measurements are more sensitive and less expensive 
than those of transitional absorption. Finally, the phase-locked four pulse wavepacket in- 
terferometry^^'^^ fit all the requirements and looks most suitable for tracking electronic and 
nuclear dynamics on the wavef unction level. This technique is applied for characterizing the 
dynamical entanglement formation in a model dimer system. 

The paper is organized as follows: Section II introduces the model of a molecular aggre- 
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gate and describes the calculation procedure. The dynamics of electronic energy transfer 
in the aggregate is discussed in Section III. Section IV presents eigenstates analysis of the 
system. Calculation of femtosecond nonlinear interferogram is described in Section V. Major 
findings of the paper are summarized in Section VI. 

II. MODEL 

A. Two-modes model of a dimer 

We consider an array (dimer, aggregate) of two coupled chromophores (monomers, 
molecules), modeled by two two-level systems being in either ground \g) or excited |e) states 
comprising for Sq and Si states of the chromophore and separated by energy ei for the first 
monomer and e[ for the second one. For convenience we choose the words "donor" and 
"acceptor" are chosen to specify the chromopfore that donates and accepts the excitation, 
respectivly. 

Since the donor and acceptor chromophores are coupled to the intramolecular vibrations 
and collective nuclear modes of the environment, the equilibrium configuration of aforemen- 
tioned modes depend on the electronic state of the dimer as illustrated in Fig. 1. Only two 
modes are taken due to following reason: According to Forster^, for a dimer, the configura- 
tion of the nuclear subsystem is characterized by two Franck-Condon active modes qa and 

with frequencies of the order of benzene stretch mode. Note that in multidimensional 
configuration space representing orientations, vibrational and other degrees of freedom we 
choose just the elongations (and stretches) that accompany the change of equilibrium; and 
refer to them as to two reaction coordinates. 

The difference in the equilibrium elongations of nuclear coordinates, corresponding to the 
ground |0) = \ga) \gb) and excited |1) = \qa) \ea) (donor), = \ga) \eb) (acceptor) states of 
the dimer characterize the strength of the electron-phonon coupling^^. For doubly excited 
state 1 2) — \ea) \eh) bath modes are elongated equally. A symbol d stands for the value of 
this elongation. The Hamiltonian of this dimer complex reads 

H = |0) Ho (0| + \l)Hi (1| + Hy (I'l + |2) H2 (2| + D, (1) 
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Here = J{|1') (1| + |1) (l'|} stands for dipole-dipole coupling, Hj for nuclear Hamiltonians: 

Here we assume that potential energy surfaces Vj{qa,qb) are harmonic and have the same 
frequency cumb and mass m for each mode and each state^°. 

^0 = [Qa + Qb) , 

^^'^ Is ri2 2\ 

^^1 = ei + (w« ~ "J +1b) ^ 

"^"^^ / 2 r ,12\ 

Vl' = €i' + (g„ + [qb-d\ ) , 

^2 = 62 + ^ (k - + h - d?) , (3) 

The expression for reorganization energy reads 

A = , (4) 

and also referred to as "Franck- Condon Energy" Epc- Figure 2 displays these potential 
surfaces Vj{qa,qb) and ground state nuclear wavepacket promoted to the state |1) by one 
ultrashort pulse from the sequence 

Vj{t) = -fiEiit), 

Eiit) = eiAi{t-ti)cos[ni{t-ti) + <^>i], (5) 

where eV, Aj, tj, flj, and stand for pulse polarization, envelope function, arrival time, 
frequency, and phase, respectively. Vj and Ej symbolize interaction energy and laser field 
strength. The dimer electronic dipole moment 

A = (|1) (0| + |2) (I'D + ^lb (0| + |2) (1|) + H.c. (6) 

allow transitions in which the exciton number changes by one. Here it is assumed that there is 
no orientational disorder and the molecular dipoles /!„ and fib are not parallel, so that pulses 
of different polarization can selectively address donor |1) or acceptor state. Restricting 
ourselves by rotating wave approximation and narrow envelope limit Aj ~ 5{t — ti), one 
can account for the first order of the laser pulse - dimer interaction resulting in the pulse 
propagation operator 

oo 



Here / labels the pulse in the sequence and equals to A, B, ... for first, second, ... pulse in a 
sequence. Subscript "x" or "y" of pulse label denotes its hnear polarization that match 
and respectivly. Note, that /7a and jib also must not be perpendicular in order to allow 
for the dipole-dipole transition J. 

B. Calculation of quantum dynamics 

After the Gaussian nuclear state has been promoted to the donor surface it starts to 
evolve in time with possibihty of the transfer to the acceptor. We calculate this dynaics 
on quantum level For the sake of convenience the dimer Hamiltonian (1) is rewritten in the 
basis of harmonic oscillators eigenstates. 



Here Mj, Nj stand for vibrational quantum numbers, \jMjNj) for nonadiabadic eigenstates 
of the dimer refered to as "diabatic", Fpc'(l'Mi/A'^i/; IMiA^^i) for Franck-Condon factors, 
describing overlaps of the wavefunctions that belong to different potential surfaces. 

The diagonahzation of the Hamiltonian (8) gives the set oik — L./cmax eigenenergies {A^} 
and the set of k^ax relevant eigenvectors combined in a form of the transfer matrix T^. 
The selected column I of this matrix gives the elements of the l-th eigenvector = (if)!^ in 
the diabatic basis. So that the solution of Schrodinger equation in diabatic state reads: 



I 

Here ip^^^{0) stands for the wavefunction at the initial moment of time expressed in the 
eigenstate basis: 



Hj = ej + nu{l/2 + Nj + 1/2 + Mj), 

D = J^^^^Fpcil'MyNir, IMiNi) llMiA^i)) {l'M[N[ 



(8) 



N^, Ml Ni 




(9) 



(10) 



Here the wavefunction is represented through a diabatic state expansion 




(11) 



where three indices can be combined into one "superindex" k. 



(12) 
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here only j — for state j — 1 for state Nj counts for number of vibrational quanta 
in a-mode, Mj counts for the number of vibrational excitations in b-mode. The maximal 
number of vibrational excitations A^max varied between 17 and 20. To achieve the highest 
numerical precision at shorter computational time we have applied so called cut-off ansatz 
to the vibrational basis set. Namely M + N < N, N -cutoff limit, M, N stand for the 
number of quanta in g^, qj, modes, respectively. 

Unless otherwise stated the diabatic initial wavefunction is taken to be coherent state in 
the donor potential well: 

Where a, (5 stand for amplitudes of coherent states in a-mode and b-mode, whose initial 
values are characterized by amplitudes |q;|, \(5\ and phases 0^. In most cases tha am- 
phtudes are expressed as integer multiples of 5 = ^J^d. Since there are four potentials, 

where one can define a coherent state, it is important to have a unified description: The 
amplitudes afi are defined so that o; = 0, /3 = corresponds to ground vibrational state in 
this potential. 

The definition of q;,/3 depends on potential surface. The mean coordinate of the 
wavepacket does not depend on potential (on electronic state). The amplitude of coherent 
state tto; /^o in ground potential has one-meaning correspondence with the mean coordinate 
qa-(lh- III order to get a unified description one may represent the coherent state of any 
potential in the basis of the ground potential; by adding the displacement between relevant 
potentials in amplitude space 5 — ^J'^d. 

The transfer matrix in Eq. 9 requirres some comments: The upper index I enumerates the 
eigenstates, while the lower index comprises for diabatic one-exciton states and combines 
three quantum numbers as defined in Eq. 12: 

C. Output variables 

The transfer of the electronic population to acceptor is 

^I'W = ^Y^i'VMNi'VMN- (14) 
MAT 
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In order to study the wavepacket interferometry signal, we have convoluted the wavef unctions 
■0, prepared by different pulses: 

CpA ^ y^^ipi'MNipi'MN- (15) 

MAT 

The donor-, acceptor-, or the whole one-exciton wavef unction are 

'^j{(ia,qb) = 01 {qa,qb\ 

^l-exciton(^^^^^) = V'l+V'l'- (17) 

Here 1-D harmonic oscillator eigenfunctions in coordinate representation for j'-th potential 
surface are denoted H'^.{qa) and H'^-'.{qb) for a- mode and b-mode, respectively. 



III. DYNAMICS 



A. Elementary act of transfer 

The simplest scenario of the electron energy transfer dynamics takes place after the dimer 
gets excited by the short pulse with narrow envelope function A{t), so that |0) -surface ground 
vibrational state is safely translated up to one-exciton surface. The pulse polarization e is 
specifically matched e _L flh to the dimer transition dipoles so that only donor surface gets 
excited into the Franck-Condon region, as shown in Fig. 3. 

Since donor surface minimum is shifted just along though the donor wavepacket 
|'0i(?o, ?&)|^ starts oscillations along this coordinate. At the time At wavepacket center 
comes closer to the "ridge" region, defined by 

e' — ei 

Qb = Qa - ^— TT' (18) 
mw^d 

The correspondent Franck-Condon window provides that part of the donor wavepacket am- 
plitude is transferred to the acceptor surface. The wavefunction amplitude in the acceptor 
state grows by small increment, linearly proportional to the intensity of the dipole-dipole 
coupling J. The transferred portion of the wavepacket, ■0i(?a, Qb) maintains the mean coor- 
dinate and momentum at the time of the elementary act, but, afterwards the motion of this 
portion of the amplitude is governed by the nuclear Hamiltonian Hi' . 
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(19) 



(20) 



performs elliptical motion about the acceptor potential surface minimum. During the time 
period that donor wavepacket stays apart from the ridge region there is no essential transfer 
of the amplitude, so the acceptor population kinetics remains fiad, parallel to the time axis. 
It can be also calculated as follows: 



B. Stepwise population dynamics 

Since nuclear potentials are harmonic, the donor wavepacket performs cyclic motion with 

period Tyib = and comes to the ridge region regularly, once per vibrational period. 

^vib 

Therefore, the elementary act of electronic energy transfer takes place repeatedly, once per 
vibrational period, as shown in Fig. 3. It is generally expected that quantum evolution of 
the coupled electronic states whose mutual detuning or coupling are modulated displays the 
stepwise character of population dynamics^^. For the short time or small coupling limit the 
almost equal portions of the wavefunction amplitude are transferred per vibrational period, 
therefore the acceptor wavefunction amplitude growth linear in time, but overall population 
growth of acceptor has a quadratic character 




a = 



{l',My,Ny\d\l',M, N) , 




(21) 



N 



Pv{t) ~ l/2(exp[ 



lA 



(22) 



for the short time limit t -C ttt. 



C. Detunings 



Fig. 5 represents the acceptor state population kinetics Pv{t) for slightly different site 
energies ei ~ ey . During first vibrational period the kinetics are indistinguishable. At longer 
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times the off-resonant energy configurations provide higher frequency of electronic nutations 
(population oscillations) and diminishes their amplitude so that acceptor population never 
gets fully populated. 

This result is in qualitative agreement with two coupled levels behavior Pda = ^ sin^ ^r^, 
where Rabi frequency JIr = (ei — e^y + 4J2 growth with detuning, and amplitude A = 
^ decreases with detuning. 

R, 

D. Population dynamics at long time limit. Revivals 

For the long time limit one expects clear and simple behavior of the population dynamics 
based on the extrapolation of the result for two coupled electronic states modeP^'^^'^^'^^, i.e. 
coherent oscillations of the population from donor to acceptor and back with Rabi frequency 

= ^/{el^^e[iy^^-^J^ . However, we show in Fig. 6 that these oscillations dephase quickly 
to the state where donor and acceptor are equally populated Pi — Py — 1/2. This quasi- 
damping originates from the destructive interference: More specifically, each single level of 
donor potential (labelled by superindex k < A^^ax) is coupled to different level of acceptor 
potential {I > N^^^). 

The coupling strength J x FC(/, k) differs for each pair. As long as many donor dia- 
batic states are initially populated, (see Eq. 13), so that the total population of acceptor 
Py — y],~^fu2 IV'/P is constructed from the sum of many contributing terms (acceptor levels 
populations) , 

2 

Sin'' {nikt)\MO)\^ (23) 

oscillating with different frequencies 

nik = ^{Hkk-Huf + AJ^¥C\\,k). (24) 

This type of dynamics was originally revealed for a two-level atom resonantly coupled to 
one-mode electromagnetic field^^. Inspite of different physical nature and different coupling 
operator the electronic population dynamics that has been calculated in this work can be 
fitted to the Jaynes-Cummings analytical formula 

Pv = e-'"'' E ^ [gVN+lt) . (25) 

N 
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\m\ 



2JFC(/,A;) 



Ik 



Here g — -^f J"^ A? / uj^nh stands for analog of Jaynes-Cummings coupling strength. The calcu- 
lated and empirical curves do coincide within collapse time interval. However, the revival of 
population difference occurs at different times for energy transfer system^^. Another differ- 
ence is that energy transfer population changes by periodic steps, as shown before, in Fig. 4. 
These studies have close association with the numerical simulation on two-mode-field JCM 
modeP^ and with the Jahn- Teller effect^®. 

For finite system, the behavior has well defined features, so there is reason to look for an 
analytical solution of exciton transfer dynamics in form 

\^l;{t)) = e-i{H,+H,,+b)t 1^(0)^ ^ (26) 
by taking in to account commutation relations between i^i, Hy, and D. 



E. Shrinking of mean trajectory 

As shown in Fig. 7 the mal transfer region Eq. 18 determines the shape of the mean 
position trajectory of the target wavepacket. Starting close to the position qa = 0,qb = 
the trajectories oscillate in both qa and q^, thought the amplitude of these oscillations in 
the "direction of transfer" growth with energy difference ei — ei/. At the time cut = n and 
cut — Stt each trajectory comes throught the same point. 

Each time donor comes to the ridge region the wavefunction portion peeled to the acceptor 
potential is not exactly the same. Each cycle acceptor wavepacket spreads wider and wider. 
This spreading makes an imprint on the acceptor wavepacket mean position trajectory, as 
shown in Fig. 7. For equal site ei — el' = configuration this trajectory repeats without 
any changes in the direction 

Q\\^ ^{Qa + Qb) (27) 
and shrinks the amplitude of the oscillations along the line 

Q± = -^(Qa - Qb) (28) 

that connects the minimas of the potential surfaces. 

The direction of trajectory shrinking depends on energy configuration, therefore it is an 
open question, whether the complete problem can be reformulated with just one vibrational 
mode, namley q±. 
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To conclude this section we have three main findings related to the dynamical behavior 
of the system: Population transfer displays stepwise character in the short time limit and 
coolapse- revival character on the long-time scale. The mean acceptor trajectory shrinks in 
amplitude along the direction depending on site energy configuration ("line of the transfer"). 

IV. STATIC FEATURES 

A. The dependence of the population transfer on the vibrational trajectory 

In our model dimer, either one or series of ultrashort polarized pulses is able to excite the 
donor potential surface into some two-dimensional coherent state Eq. 13. For example: As 
shown before, a single x-polarized pulse creates the coherent excitation in the Qa mode having 
Pranck-Condon amount of vibrational energy, placed initially at [qa = 0,qb = 0) with initial 
phase (f)a = 0. However, the specific sequences of pulses can create vibrational excitations in 
donor surface, starting at different points of phase space. One of these examples is illustated 
below: Let's consider the excitation prepared by a y-polarized pulse. After a quarter of 
vibrational period Tvib/4 we send an x-pulse and, another quarter period later, a y-pulse. . 
This pulse sequence creates a two-dimensional circular motion about center of donor surface, 
having AEpc energy in both modes moving 2d apart from the center of donor surface, starting 
at point Qa = d, qh = 2d. When the last pulse is applied half a period (rvib/2) later then 
it produces again the qa coherent excitation in a donor surface having, however no initial 
elongation but nonzero momentum giving an initial phase 0^ = 37r/2. This section shows 
that specific pulse series does create specific coherent states in a donor surface. 

The change of vibrational states in donor surface does affect the intensity of the electronic 
energy transfer to acceptor. To investigate this we have considered a set of coherent states 
having the same amount of vibrational energy {Epc) differently apportioned between qa 
and qb modes, see Fig. 8. As shown before, the elementary act of transfer takes place when 
the wavepacket crosses the ridge region of potential energy landscape Eq. 18. The mean 
coordinate trajectories of different coherent states cross this line in a different manner. In 
accordance with Landau-Zener formula^^'^^, as longer the wavepacket stays on the ridge fine 
as faster the population transfer goes. 

The minimal intensity of the transfer is found for the intitial state of donor having no 



11 



vibrational excitation at all (ground vibrational state). This initial state provides simple 
oscillations of electronic amplitude from donor |1) to acceptor and back with frequency 
J X FC(1, 0, 0; 1', 0, 0) in leading order determined by Pranck-Condon overlap of ground 
vibrational state of each potential. Vibrational trajectory in q\\ direction provides similar 
oscillatory behavior on the long time scale. Excitations of q^,., q±, and circular excitation 
lead to the quicker transfer on short time-scale and collapse-revivals behavior on the long 
time scale. It is clearly shown that presence of vibrational excitation enhances the transfer 
of electronic amplitude between one-exciton states. 

B. Origins of parallel, perpendicular and combined effect 

The ridge region is intersecting the line connecting the minima of donor (g^ = d, qh = 0) 
and acceptor (g^ = 0, qi, = d) surfaces. It is reasonable to measure the distance between 
wavepacket and the region of the optimal transfer along this line. Therefore we use the 
45 degrees rotated system of coordinates consisting of q\\ and g_L coordinates, described by 
Eq. 27 and Eq. 28. The motion of a coherent wavepacket along q^^ is expected to determine 
the efficiency of the transfer. The dependence of population transfer on the motion in g|| 
does remain to be an open question. Instead of taking an exhaustive collection of all possible 
two-dimensional states, in Eq. 13 we consider the set of vibrational states having different 
amount of excitation and phase along either or gy, in order to exploit all available transfer 
regimes. 

C. Dependence on energy difference 

As long as region of potentials' intersection location Eq. 18 qb = qa — ^^J^^d depends on 
site-energy difference between donor and acceptor moieties, the efficiency of the population 
transfer is expected to depend on the difference ei — — E^a- The dependence on energy 
difference gives a sence how The energy Specific rolecules in specific solvents give various 
regimes of energy difference. By scanning all values of energy difference we get a sence of 
behavior of various real systems. 
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1. The role of pendicular excitation. 



It follows from Fig. 8 that the presence of vibrational energy in q± - mode gives rise to the 
population transfer. There is also no strong dependence on energy difference E^a — e^ — 
because almost any position of the intersection line Eq. 18 Qh — % — (ci — eir)/mu!^d is 
reachable by g^-coherent wavepacket. The intersting dependence of population transfer on 
the phase of q± coherent motion is left to consider later. 

D. Dependence on parallel excitation 

1. Introduction of Fock states 

The dependence of population transfer on the excitation of the g|| mode is rather small. 
Therefore,we do not consider the initial phase of of the q\\ excitation but only the amount of 
vibrational energy in this mode. The relevant state with the definite number of vibrational 

quanta is so-called Fock-state with P vibrational quanta in q\\ mode. In the basis of the 
natural vibrational quantum numbers M, N for the modes Qa, Qb this state reads: 



Here {qa\ |0) = '00 (?a) {Qb\ |0) = i^oiQb) stand for ground states in each vibrational modes. 
Since the common ground state is factor of those two, one gets two dimensional M, A^- 
Fock states in q\\, q± by applying the derivatives along these coordinates M and N times, 
respectively. 



E. Overall Markus' hump analysis 

The simplest vibrational state is the ground state, having no vibrational quanta at all. For 
this state the dependence of acceptor population on energy difference E^a — ^i — has no 
admixture of vibrational influence, as shown in Fig. 9. This dependence has a form of overall 
hump, modulated by fringe-like srtucture. The overall hump has maximum at Eda — ^i — 




(29) 




(30) 
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ei' — —2Efc, where the acceptor diabatic potential surface crosses the minimum of the donor 
potential surface. This corresponds to the activationless regime of the electron transfer with 
one reaction coordinate in the Marcus theory, which has an enormous range of applications to 
exciton, electron, proton transfer and many other chemical reactions^'^^'^^'^^'^^. The fringes 
originate from the individual resonances between vibrational levels, belonging to donor and 
acceptor moieties. The presence of theese individual resonances support the discussion in 
section III C and Eq. 23. Up to our knowlege such fringes were at first noted for one-mode 
system^^. 

For large value of energy difference, vibrational Fock-state in the gy-mode provides faster 
population transfer than the ground vibrational state of the donor potential. In order to 
understand this effect an cigcnstatc analysis has been performed. As far as one knows 
the vibronic eigenstates of the dimer, it is possible to find the mean values of some relevant 
variables, (like e.g. mean coordinate in Eq. 21). For example, the mean values of momentum 
P|| = m^ii and q± — 'mq± are calculated systematically for all eigenstates and develop a 
regular structure displayed in Fig. 10. 

The diabatic ground state and a Fock-state in donor surface were expanded over eigen- 
states basis set 

|^Sround>= J]Pr"">>, (31) 

|V'Focfe>= J]Pj'°^1^>- (32) 

and displayed in Fig. 10. Here \u) stands for u-ih eigenstate. The set of mean momenta 
enumerates eigenstates. The ground diabatic state involves the eigenstate that has minimal 
mean momenta. The "parallel Fock state" does not have any vibrational excitations in q± 
direction. That is why it is expected to employ just those eigenstates with larger mean values 
of In contrast, the numerical simulation shows, that this state employ some eigenstates 
with large momentum p± in "perpendicular" direction. The presence of such states in 
the eigenstate expansion of the Fock-state is, probably, responsible for the difference in 
population transfer rates, provided by these two diabatic vibrational states. 
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F. The role of perpendicular phase 



1. Specific values of phase 

We return back to the dependence on excitation of q± mode. The challenging question 
is whether the transfer rate depends on the amount of vibrational energy in this mode only 
or not. Alternativly, it can depend on the initial phase of the q± -coherent excitation in 
donor manifold. We demonstrate the results for two most distinctive cases: Wavepacket is 
far apart from the intersection ridge = 0) and the opposite position (0_|_ = tt) that stays 
at closest to the acceptor surface minimum, a — Se^^ , (3 = de~^'^ and ... 

2. Evident results 

Figure 11 shows that there is small but evident distinction in population transfer, corre- 
sponding to this two cases. For positive site-energy difference Ejja — — > referred 
to as "uphill transfer" the coherent state initially positioned in the acceptor region provides 
quicker transport. For negative site energy difference difference Eda = t\ — ey < 0, referred 
to as "downhill transfer" the quicker transfer occurs to be for the coherent state, initially 
placed in the region of donor surface minimum. 

G. Explanation to repositioning of adiabatic potential 

The difference in the transfer rate becomes transparent while analyzing wavepackets' 
motion in the adiabatic potentials, as shown in Fig. 12. The analysis is based on the fact 
that the actual wavepacket motion takes place in the upper and lower adiabatic surfaces 



Therefore, the time evolution of this system depends on how the initial wavepacket is rcpar- 
titioned between upper z/+ and lower u_ adiabatic potentials. For given mean position and 
mean energy of wavepacket the largest part belong to the surface that have the closer energy 
value at this position. For any energy difference the wavepacket (f)± — belongs to the lower 
potential, while the wavepacket with shifted phase 0_l = tt belongs to the upper potential 
and is confined there. 




(33) 
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The only difference between "uphill" and "downhill" transfer is the position of the min- 
imum of the upper adiabatic potential. As shown in Fig. 12 this minimum is shifted aside 
of donor or acceptor, for downhill and uphill cases, respectively. 

For downhill transfer the 0_l = vr wavcpacket stays longer in donor region, therefore the 
transfer to acceptor is diminished. For uphill transfer the 0_l = vr wavepacket stays longer in 
acceptor region and so provides the faster transfer with respect to the transfer corresponding 
to the (f)± — coherent vibration. 

In this section has been revealed that energy transfer depends on regime of launched 
wavepacket passing through the cotential-crossing region of energy landscape mostly asso- 
ciated with amplitude of motion along the line connecting the potential minima. 

V. DETECTION 

A. General principles 

There has been discussed nontrivial features of energy transfer channel of photoexcitation 
(evolution / decay) in dimers that would be interesting and relevant to reflect by means 
ultrafast spectroscopical measurements. 

In this section, we discuss the novel measurement scenario, allowing to get wavefunction- 
amplitude- level information about exciton transfer and nuclear dynamics in molecular 
dimers. This scenarion requires four polarized-pulsc excitations. In order to be concrete we 
choose the specific number of pulses and polarization arrangements: four pulses AyBxCxD^ 
separated by preparation tp — ts — tA, waiting — tc — tsi and delay — to — tc times, 
respectivly^^'^^. As displayed in Figure 13, our "observable" is the /|| component of the 
fluorescence from the oriented sample that has been excited by aforeshown pulse sequence. 
The /|| fluorescence intensity is proportional to the population of the acceptor state of 
the dimer. There arc few elementary laser-molecule interaction processes that give rise to 
the population of the state. All but nonlinear (depending on all four pulses) processes can 
be cutted off out by applying mechanical choppers^^'^^. 

Among the rest of the nonlinear processes we consider only those quadrrilinear in the 
intensities of each pulse because they are of larger value 

Pv{AyB,C^D^) = 25R{ (l'| {Ay\ \JD,C,B,) e-'^--'^'^ 
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{Ay\\D,C,JB,) e-'^-+'^'^ 
{DxBxA.y\ \ JCx) 

{CxBxAy\\JDx)\l') }. (34) 

Note that in this polarisation scheme all contributions to Py — {il^v \ depend on exciton 
transfer. Linearly combining some quadrilinear fluorescence signals taken with different 
relative phase-locking angles 0^ between AyB^ and CxD^ pulses as described in^^ 

{{Ay\, \{DxCxJBx\,f^ = 
^{P.,,O.O)-P.(|.-|)+.A,(|.O)-,P.(O.0} 

one gets the value that can be taken in to account by the double side Feynman diagram^^'^^ 
shown in Fig. 14. Here the left side describes so-called "target" wavefunction, where the 
population amplitude in state |1') is promoted from the ground state by three pulses B^, 
Cx, Dx and one primariy electronic energy transfer act between Bx and Cx pulses. Here the 
vertices of the fiagram are given by linear terms of pulse propagators 7. 

The actual population of state |1') (and fluorescence from it) is proportional to the 
coincidence between "target' and so-called "reference" wavepacket, schematically presented 
by the right part of the Feynmann diagram and promoted from the ground state to the 1 1') 
by a single Ay pulse. 

As long as electronic transitions in this model dimer are coupled to the nuclear modes, the 
probability amplitude wavepackets change their positions and may not coincide for the left 
and right part of the diagram. This argument gives reason to expect essential interference 
population of |1') for some speciflc values of time delays only. 

B. Equal site energies versus downhill 

Figures 15, 16 show the results of the numerical simulation of nonlinear interferometry for 
two simple cases: equal site energies of donor and acceptor ei — ei' = and downhill transfer 
ei — Ey — 2Epc- As long as many other processes have been excluded from this signal so it 
shows negative values. However, the actual fluorescence intensity is always positive so the 
negative signal shown in this flgure means nothing but minimum of of actual fluorescence 
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intensity. The interferograms differ in three features: (i) intensity, (ii) time of maximal 
intensity, (iii) interference fringe structure. 

C. Matching conditions 

The theoretical estimation of the "preferred" delay times is based on the quasiclassical 
analysis. The mean position and momenta of target and reference wavepacket must coincide 
at time moment just after pulse. In order to find coincidence criteria we perform few 
assumptions: (1) nuclear wavepackets evolve under Hj, (2) For the sake of convinience Cx, 
Dx pulses and free evolution between them are transferred to the "reference" side of the 
convolution: 

< Ay\DxC^B^ >^< AyD+C+\B, > . (35) 

Mean while, between the pulse arrival the dimer state performs free unperturbed evolution 
(vertical solid arrows on Feynmann diagram 14) symbolized by square brackets and evolution 
time [t] — exp {—iHt) so the term of our interest contributing to interference population 
reads 




1^1' > 



Here are exact definitions of target 

\av) = e^"-*--^^''*'^ (I'l Ay[tp + + ta]D+[-ta]Ct |0) (37) 

and reference 

lei') = (l'|[t^]5,|0) (38) 

wavepackets. Note that because of bra-(l'| Eqs. 37-38 represent the nitc/ear wavefunctions in 
elecronic state 1'. While preparing reference side of interferometry population corresponding 
to process 36 the wavefunction always have the following form 

(V'l ~ 01 {a\ (39) 

Here (j| stands for electronic state 0, 1, 1', 2, {a\ and (/?| stand for coherent states Eq. 13 in 
modes and respectivly. The amplitudes 

[mw ( .Pa\ 

„= /^k + ,iL}, (40) 
V 2 I muj ) 
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determine the mean position and momentum of the coherent state. Depending on the 
electronic state {j\ the mean values a{t), (5{t) circle around the minimum of the relevant 
harmonic potential surface: (0,0) for |0), (5,0) for (0,5) for and {5,5) for |2). Since 
ip) tw, td are chosen one may follow the phase-space trajectories of a, j3 for target and ref- 
erence states as shown in Fig. 17. In contrast to the reference wavefunction being always 
in either one of electronic states the target wavepacket reflects the dynamical entanglement 
formation between |1 >, |1' >, and nuclear modes. We specialize tyj equal to half of vibra- 
tional period Tyib = 27r/u' to prevent multiple acts of the transfer. Target wavepacket evolves 
during the ty, as discussed in section III A. Initially it is nothing but coherent excitation in 
donor surface starting in position (g^ = 0, = 0) and oscillating about (g^ = 0, ^5 = rf). 

Here wc refer to the new variable A so that after At^ time interval after pulse the 
first elementary act of transfer takes place. The variable A accounts for this moment of time 
only approximately and depend on energy arrangement of the potentials Ed a — ^1 — ^v- 
The actual transfer is not instantaneous but continuos. After the act the transferred portion 
of the wavepacket amplitude does maintain its mean position and momentum, but starts to 
evolve about the minimum of the acceptor parabola (amplitudes in the basis of ground state 
potential). After it evolves in such a manner during the (1 — A)ty^ time its position and 
momentum need to coincide with mean position and momentum of the reference wavepacket. 
The analytical condition of such coincidence reads 

a - mode : (e^'""^*"' - l)e-^'"*» - 1 = -e^^*^ 

h - mode : e-^"'([i--4]*-+*'*) = ^-''^itp+tu.+ta) _ (41) 

As far as tu, is strictly fixed to be Tvib/2 the conditions Eqs. 41 can be considered as system 
of two algebraic equations in respect to two variables tp and td- The solution of this equation 
gives the estimation when the target and reference wavepackets do overlap at best: 

tp = {m + 1 - A/2) Tvib, 

td = (n + ^/2)Tvib. (42) 

It is integer multiples of vibrational period shifted far off resonance in the direction of smaller 
tp and larger td- ^ ~ is expacted to be a bit larger than zero for equal site energies and 
^ ~ 1/4 for the activationless downhill transfer. This analythis goes along with the results 
of the numerical simulations. 
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The calculated phase-space trajectories for target state are shown in Fig. 18. Starting and 
end points do not stay on the quasiclassical trajectory (on the FC energy shell) . However, 
increase of site-energy difference shifts the trajectory end further (closer) from (0, 0) point 
in a- mode (b-mode). 

D. Fringes 

Figures 15, 16 also display fringes in target nuclear wavepackets. The maximal coinci- 
dence of target nuclear wavepacket with referrence wavepacket is one probable reason of 
the interferogram fringe structure. As long as interference structure is arranged to have 
no fringes along the line tp = td but change the sign along the line tp = —td the following 
speculations can be suitable: Change of the net signal phase of the phase-locking factor 
exp —iluJiockiip ~ td) takes place along the line tp — td- This factor is zero for tp = td (and 
constant for tp — td + x). In case there is an inhomogeneity between different dimers that 
form the sample, this factor will still have maximal value for tp — td- This situation is 
analogeous to the stimulated photon echo scenario where the last pulse just plays the role 
of probe in order to detect the created net growth of the transition dipole moment mean 
value (homodync detection). 

Fringes in signals reflect (i) matching / mismatching of phase (ii) frequency of signal 
oscillations along tp and td Information from fringe structure: Fringes Frequency corre- 
sponds to the mismatch of energy configuration and locking frequency. The difference of 
continuous and instantaneous transfer. An alternative interpretation rests on the velocity 
of the vavepacket. 



TABLE 1. Position of signal amplitude maximum and rate of its phase change taken at 
the point of maximal signal amplitude. 



configuration 


_tp_ jd_ ijp r\j_ 

Tvib Tvib 


equal ei — ei' — 0, quantum 
equal ei — — 0, semiclassical 
downhill ei — ei' — 2£'fc, quantum 
downhill downhill ei — — 2Efc, semiclassical 


0.923 1.077 -0.53 + iO.Ol 0.53 - iO.Ol 
1.000 1.000 
0.750 1.250 -7.52 + iO.Ol 0.13 + iO.Ol 
0.750 1.250 -7.39 
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VI. CONCLUSIONS 



We investigated the dynamical entanglement formation in a simple dimer with two nuclear 
modes energy transfer model system and its reflection by means of ultrafast nonlinear phase- 
locked wavepacket interferometry. We have revealed some intriguing features that may 
attract an attention of physical chemistry and quantum optics communities. 

Following the ultrafast excitation of donor the population of acceptor state gets riscd by 
small increments once per vibrational period. The long-time population oscillations between 
donor and acceptor states display collapses and revivals similar to those in Jaynes-Cummings 
model. The mean coordinate of the acceptor wavepacket loses its amplitude with time. 

The intensity of the transfer is found to grow with depth of donor wavepacket penetration 
into the acceptor region. This depth depends on site-energy difference of donor and acceptor 
and amount of vibrational energy in donor potential in the transfer direction (perpendicular 
mode). The minor influence on the transfer intensity comes from initial phase of donor 
vibrational coherence and amount of vibrational energy in the parallel mode. 

The consideration of the four-pulse phase-locked nonlinear wavepacket interferometry ex- 
periment on such model dimer shows that dimers with different site-energy difference provide 
different nonlinear optical response. The nonlinear wavepacket interferograms corresponding 
to equal-site energy dimer and activationlcss dimer have maxima at different delays between 
excitation pulses. This difference is predicted with satisfactory precision using quasiclassical 
analysis of the wavepacket mean trajectories. 

The possible future development of this research should include the consideration of dis- 
oriented sample. Further work will account for vibrational relaxation and inhomogeneity of 
site energies and dipole-dipole coupling induced by various spatial orientations and distances 
between monomers, as well as application of these findings to the real molecular aggregate. 
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doubly excited |2> = leg>leg> 
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acceptor state |1'>= leg>lg|^ 
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donor state |1> = l^l^ 
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FIG. 1: The reasons to use two modes; coupled chromophores (circles) can be both in ground 
(lower plot) state, or in a one-exciton state when either one of the chromophores is excited (left 
and right plots), or doubly excited state (upper plot); The monomer's static dipoles Da, Db (thick 
arrows) change for excited states; The equilibrium orientation of inter- and intra- molecular nuclear 
configuration (e.g. solvent librations - small arrows) reorganizes in order to compensate for the 
change of Da, Db- 
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FIG. 2: Potential surfaces of dimer illustrating equation 3 (upper plot); an example of mutual 
spatial orientation between laser pulse Eq. 5 and transition dipoles of dimer's species Eq. 6 (lower 
plot) 
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0.5 1 1.5 2 

FIG. 3: Elementary act of the electronic energy transfer; Contour lines of probability density 
wavepackets: iV'iP - in donor potential (thick dashes), and IV'P - in acceptor potential (thick solid 
lines) are displayed at different time slices t = 0, t = rvib/4, and t = 

taMvib/2; inset shows the correspondent population kinetics, stars correspond to the time slices; 
ridge line qa = qt connects "0" and "2"; ei — ei' = 0, J = oj/lOO 



26 




FIG. 4: Short-time population dynamics for coupling strength J = w/lO and various energy 
configurations: ei — ey = (solid), ei — = io/2 (dashed), ei — ey = 7io (dotted), ei — = 
7.39a; = 2£'fc (dot-dashed); solid and dashed displays one step per vibrational period, dotted and 
dot-dashed dispaly two steps per vibrational period; since short time-population is small log-scale 
is used for Pi/. 
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2 4 6 8 10 12 14 16 18 



FIG. 5: Population kinetics for different off-resonance detuning: ei — ei/ = (solid) ei — ey = uj/16 
(dashed) ei — ei' = uj/8 (dot-dashed) ei — ei' = uj/4 (dots) 
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FIG. 6: long-time population dynamics for ei — ei/ = 0, J = uj^i^^iq. collapses and revivals, 
simulation (solid) and equation 25 with two sets of parameters: g = JE-pc\/^/oj, a = ■\/muj/2d 
(dashed) and g = JE^c \/2/a;, a = ^mud (dot-dashed), 
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FIG. 8: Dependence of transfer on the specifics of the vibrational excitation in donor surface: mean 
position trajectories of donor surface wavepacket (upper panel) and correspondent populations 
kinetics at short time (lower panel) and long time (inset of lower panel); ei — ei'=0;J = u/lO; for 
each trajectory the vibrational energy of Epc is differently apportioned between nuclear modes: 
Franck-Condon excitation, ga-mode excited (solid line), g_L-mode excited (dashed line,) g||-mode 
excited (dot-dashed line), circular excitation - both modes equally excited (dotted line); for the 
purpose of reference thin solid line on bottom panel correspond to ground nuclear state in donor 
potential; for better illustration qg_, q±, trajectories are prepared to have a small portion of 
vibrational energy in the adjacent coordinate. 
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Influence of Parallel Coherence 




-15 -10 -5 

FIG. 9: Electronic energy transfer to acceptor state for various site energy differences for J = oo/2] 
solid line with circles stands for zero-point vibrations in donor surface (no vibrations), dashed line 

- Fock state excitation in donor surface, upper panel shows population of at t = 2rvib, lower 

_ °° 

panel shows acceptor population averaged over T = lOOTvib: Pi' = ^ / Pv{t)dt 
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FIG. 10: Mean values of momenta for 1-exciton eigenstates \l); each dot correspond to one 
eigenstate; horizontal position corresponds to the value of |/); vertical position corresponds 

to the value of {l\p\\ \ \l) here ei — ey = 0, J = cj/lO; The J-coupling strength induces splitting in 
the direction of transfer. 
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cut-off ansatz, Q^^^^ coherency, J=w Ep =3.69co^^,, N^,,=20 




FIG. 11: Electronic energy transfer to acceptor state |1') as function of energy configuration 
calculated for initial vibrational excitation of g_L-mode having Epc of vibrational energy in donor 
potential; initially posiitoned apart from acceptor ((p = 0, solid), or close to acceptor (</> = vr, 
dashes); for ei — ei' < = provides faster transfer than </> = vr for ei — ei' < = provides 
slower transfer than cp = tt. 
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FIG. 12: Schematic illustration of adiabatic potentialenergy surfaces Eq. 33 for ei — ei' < (left 
panel) and for ei — ei' < (right panel) solid line stands for lower potential dashed line stands 
for upper panel i^^; gaussian profile symbolizes the initial vibrational state in donor potential 
surface: having zero phase of q± excitation (left), or having vr phase of q± excitation (right). 
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FIG. 13: Schematic representation of interferometry experiment; ensamble of dimers with specific 
spatial orientation (of transition dipoles /Udonor and /Xacceptor) is labelled as "sample"; "laser setup" 
generates four short pulses, resonant to dimer excitation frequency, pulses are labelled with capital 
letters, subscript "x" or "y" indicates pulse polarization; time delays between pulses are labelled as 
tp, iw) t(i, where subscripts p, w, d abbreviate for "preparation", "waiting", and "delay", respec- 
tivly; excited sample's polarized fluorescence is measured with "detector"; only the fluorescence 
part with polarization attributed to "acceptor" chromofore is analyzed in the text. 
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FIG. 14: Double-side Feynmann diagram, reprsenting processes contributing to interference pop- 
ulation Py =< 0\Ay\ X |DxCx-Bx|0 >; left part symbolizes ket-vector \DxCxBx\0 >; right part 
symbolizes bra- vector < 0|^y|; straight arrowed line stands for time evolution of a single chro- 
mophore (two- level-system); wiggly line stands for acquanted (tail down) or emitted (tail up) 
photon; doubly wiggly line stands for exciton transfer between chromophores. 



37 



1.5 

equal site, 




p 



FIG. 15: Real part of the interferometry signal provided by {AyD^Cxl \ JBx) - process for ei — ei' = 
energy configuration as function of delay times tp, td between pulses; dark areas correspond to 
negative signal, bright areas correspond to positive signal, contour dashed lines display the absolute 
value of the signa; the signal has main maximum aleft far of the diagonal repeating periodically at 
integer multiples of tp, with slow change of the sign yelding three fringes: oscillation frequency 
—0.53a; along tp, 0.53ia; along td] there is also small satellite peak at tp = 0.5rvib, td = 0.5rvib- 
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FIG. 16: Real part ofthe interferometry signal provided by {AyDxCx\ \ JBx) - process for ei — 
ei' = 2Eyc enegry configuration as function of delay times tp, td between bulses A, B and C, D, 
respectivly; dark - negative signal, bright - positive signal, dashed contours - amplitude of signal; 
the signal maximum at tp = 0.75r, = 1.25r repeats periodically with slow change of phase along 
t^-axis: — O.SSo;, and relativly fast phase change along tp-axis: +7.53w; There are also two satellite 
peaks at tp = 0.5t, = 0.5r and tp = t, t^ = t. 
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FIG. 17: Scheme of quasiclassical matching, a-mode and b-mode phase-space wavepacket center 
trajectories for acceptor potential; Left panels illustrate reference wavepacket Eq. 38, a-mode: 
a = (5 (l -I- e~*'^*^*p+*"')) (upper panel), b-mode: /? = 5 (l -|- e*'^*'*) (lowe panel), dashed circles 
symbolize Pranck-Condon Energy shell; Right panels illustrate target wavepacket Eq. 37 a-mode: 
a = 6 {e'^^'^-A — 1^ g-«^<™ (upper panel), b-mode: (i = 5 {l — e^^i'^-A-tw)'^ (lowe panel); two arrows 
on lower panel symbolize the wavepacket center motion before (first arrow) and after (second arrow) 
the instant of transfer. 
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FIG. 18: Calculated phase space trajectories for target wavepacket for ei — er = (solid), ei — ey = 

Epc (dashed), ei — ei/ = 2Epc (dotted); during the time interval < tot < tt, dot-dashed line points 
quasiclassical phase space trajectory in dispalced harmonic potential (Franck-Condon energy shell) 



41 



1 y I I I I I 




